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ABSTRACT 

We examine the connections between the current global star formation activity, luminosity, dynam- 
ical mass and morphology of galaxies in the Local Volume, using Ha data from the 11 Mpc Ha and 
Ultraviolet Galaxy Survey (IIHUGS). Taking the equivalent width (EW) of the Ha emission line as 
a tracer of the specific star formation rate, we analyze the distribution of galaxies in the M^-EW 
and rotational velocity (Kna£c)-EW planes. Star-forming galaxies show two characteristic transitions 
in these planes. A narrowing of the galaxy locus occurs at Mb '-^ — 15 and Vmax ^ 50 km s~^, where 
the scatter in the logarithmic EWs drops by a factor of two as the luminosities/masses increase, and 
galaxy morphologies shift from predominately irregular to late-type spiral. Another transition occurs 
at Mb ~ —19 and Vmax ~120 km s~^, above which the sequence turns off toward lower EWs and 
becomes mostly populated by intermediate and early- type bulge-prominent spirals. Between these 
two transitions, the mean logarithmic EW appears to remain constant at 30A. We comment on how 
these features reflect established empirical relationships, and provide clues for identifying the large- 
scale physical processes that both drive and regulate star formation, with emphasis on the low-mass 
galaxies that dominate our approximately volume-limited sample. 

Subject headings: galaxies: dwarf- galaxies: evolution - galaxies: fundamental parameters - galaxies: 
ISM - galaxies: stellar content - stars: formation 



1. INTRODUCTION 

It is well established that many of the global prop- 
erties of present-day galaxies scale coherently with one 
another. Of particular interest is the star formation ac- 
tivity; understanding and deciphering trends with this 
property provides insight into the astrophysics that con- 
trol the pace and extent of a galaxy's evolution. Perhaps 
the most well-known trend is that the current star for- 
mation rate (SFR) varies monotonically along the Hub- 
ble sequence (Roberts 1963; Kennicutt 1998). The scale 
begins at the ellipticals, which generally have red col- 
ors, little gas and essentially no star formation, and 
proceeds through the lenticulars, spirals and irregulars, 
which progressively have bluer colors, larger gas fractions 
and greater relative levels of star formation. 

The quantification of these patterns began with the de- 
velopment of SFR measurement diagnostics (e.g., Tinsley 
1968; Kennicutt 1983a), and is now reaching full matu- 
rity with the application of such techniques to large, sta- 
tistically complete samples of local galaxies (e.g., Brinch- 
mann et al. 2004, Salim et al. 2007), as well as of galaxies 
at higher redshift (e.g. Noekse et al. 2007a). Progress 
has been rapid in quantifying relations for galaxies in 
the intermediate and high-end of the luminosity func- 
tion, with the significant result that the galaxy distribu- 
tion is bimodal and may divide into distinct populations 
of old, red spheroids and blue star forming disks at a 
stellar mass of - 3 x IO^^Mq (Kauffmann et al. 2003). 
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In contrast, progress on the dwarf galaxy population 
has been slower, primarily due to the lack of complete 
and statistically robust samples in the low-mass regime. 
To address this need, we have focused on our near- 
est neighbors and are taking an appoximately volume- 
limited (and hence dwarf-dominated) inventory of star 
formation in the Local Volume. Our 11 Mpc Ha and 
UV Galaxy Survey (IIHUGS) has obtained narrowband 
Ha emission-line imaging, which traces massive 0-star 
formation, and is also collecting GALEX UV imaging, 
which traces the photospheric emission of 0-type as well 
as longer-lived B-stars. IIHUGS also serves as the basis 
for the Spitzer Local Volume Legacy program, an IRAC 
(mid-IR) and MIPS (far-IR) imaging survey of 258 galax- 
ies, which will capture the energy which is absorbed and 
re-radiated by dust. The IIHUGS parent sample con- 
tains ^^400 galaxies and includes all known spirals, irreg- 
ulars, as well as star- forming early- type galaxies with d < 
11 Mpc, \b\ > 20 °and B < 15 mag. This Letter takes 
a first look at results from the recently completed Ha 
imaging component of the survey (Lee 2006, Kennicutt 
et al. 2008), and explores trends in the magnitude vs. 
Ha equivalent width (EW) plane, a close cousin of the 
more commonly studied color-magnitude diagram (e.g. 
Baldry et al. 2004 and references therein). 

2. DATA 

To visualize trends in star formation across the 
IIHUGS sample, we examine the relationship of the Ha 
EW with two indicators of the mass, Vmax , the maximum 
circular rotational velocity, and Mb- We primarily use 
Mb since i3-band photometry is widely available in the 
literature and can be approximately reduced to a com- 
mon (RC3) system for all of the galaxies in the sample. 
A well-known drawback with using the i?-band light to 
trace stellar mass is that the mass-to-light ratio in the 
blue has a strong systematic dependence on the recent 
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We provide details on the IIHUGS sample selec- 
tion, observations, data processing, calibrations, mea- 
surements and completeness properties in Lee (2006), 
Kennicutt et al. (2008) and Lee & Kennicutt (2008). 



Fig. 1. — The Local Volume star-forming galaxy sequence in the (a) Mg-EW and the (b) rotational velocity-EW planes. Galaxies in 
the core sample of IIHUGS (i.e. those with T>0, D<11 Mpc and \b\ > 20°) are shown. Gray dashed lines are drawn at AIb = —19 and 
V'max=120 km s~^ , and Mb = —15 and Vmax=50 km s^^ to indicate the two transition regions discussed in the text. 

star formation history itself (e.g. Bell & deJong 2001). 
Examining the light in redder passbands, particularly in 
the near- and mid-IR, mitigates against these effects, and 
the 3.6 /im and 4.5 /im imaging being collected by the 
Local Volume Legacy will especially useful in this regard. 
In the meantime, we have also looked at the variation of 
the EW with Vmax as a qualitative consistency check. 
Vmax traces mass gravitationally, and should in principle 
be independent of the recent star formation activity in 
individual galaxies. The existence of the baryonic TuUy- 
Fisher relation over the full mass-spectrum of late-type 
galaxies (e.g., McGaugh 2005) supports the use of Vmax 
as a mass indicator even in our dwarf dominated sample. 
Using Vmax is not without its own weaknesses however. 



One issue is that the maximal rotation velocity for low- 
mass systems is somewhat ill-defined since dwarfs are 
generally known to be solid body rotators (e.g. Skillman 
1996). Further, previous measurements of Vmax arc lim- 
ited, and are available for only ^75% of the sample. Here 
we use Vmax values taken from the HyperLeda database, 
which have been extrapolated from single-dish 21-cm line 
widths (Paturel et al. 2003) using galaxies which have 
both line widths and rotation curves in the literature. 

To track the normalized star formation activity, we use 
the Ha EWs from IIHUGS. Recall that Balmer emission 
results from the recombination of nebular gas which has 
been photoionized by young, short-lived O and B stars. 
Thus, the EW, which is given by the Ha flux divided 
by the red continuum flux density, is an indicator of the 
strength of the current SFR relative to the total mass 
of stars (i.e., the specific SFR), and of a closely related 
quantity, the ratio of the current SFR to the past av- 
erage SFR (i.e., the stellar birthrate parameter b). The 
IIHUGS EWs used here have all been measured from 
narrowband imaging and are integrated over the entire 
extent of the galaxies so they are representative of global, 
galaxy- wide averaged values. They do not suffer from the 
aperture effects that plague fiber-based spectral measure- 
ments. 



3. RESULTS 

The distribution of IIHUGS galaxies in the Mb-EW 
and T/max-EW planes are shown in Fig. 1. Mb is based 
on distances calculated from flow-corrected velocities as- 
suming Ho=75 km s~^ Mpc~^ when no direct distance 
estimates are available. No internal extinction correc- 
tions have been made to Mb- We plot only those ^300 
galaxies in the sample that are (i) on the late Hubble se- 
quence beginning with the SOa morphological type (T>0) 
and (ii) outside the Galactic plane (|&| > 20°), and (iii) 
within 11 Mpc. It is within these limits that we have 
tried to be as complete as possible in our compilation of 
currently known galaxies. This core sample of IIHUGS 
is statistically complete down to B=15 (Lee 2006, Lee 
& Kennicutt 2008). Although no corrections for internal 
extinction and the [Nil] contribution to the EW have 
been applied to the data shown here, we have checked 
that such corrections (as described in Lee 2006 and Ken- 
nicutt et al. 2008) do not affect the results that fol- 
low. We have used different symbols to represent differ- 
ent morphological types (as given in the RG3) as indi- 
cated in the plot. Upper-limit symbols represent galaxies 
undetected in Ha. In Fig. lb, galaxies that are nearly 
face-on (i < 30°) are excluded because large inclination 
corrections result in highly uncertain Vmax values. 

Over the full range of luminosities and circular veloci- 
ties, star forming galaxies trace a continuous sequence in 
morphological type and EW with two characteristic tran- 
sitions. The trends are qualitatively consistent in both 
diagrams which gives us confidence that the features are 
real and not merely due to covariance between Mb and 
the EW, and/or biases in the available Vmax data or its 
measurement. 

Dwarf galaxies, mostly Magellanic spirals and irregu- 
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lars, appear as the broad swath of light blue symbols at 
the bottom of the diagram. They span the largest range 
of EWs observed in galaxies, and include both high EW 
starbursting blue compact dwarfs, as well as low EW 
systems that are barely forming stars and thought to be 
in an intermediate phase between the the gas-rich dwarf 
irregulars and gas-poor dwarf spheroidals. In contrast, 
above ~ 50 km s~^ and Mb ~ —15, a transition ap- 
pears to occur where the EW distribution tightens. This 
change is accompanied by a shift in the dominance of ir- 
regular morphologies to late-type spirals in the sequence. 
Galaxies in the "waist" of the sequence have an average 
EW of ~30 A. Moreover, galaxies in which star forma- 
tion has turned off or which are undergoing extreme star- 
bursts are virtually absent. 

Finally, there appears to be another transition at ^ 120 
km s^^ and Mb ^ —19. Above this velocity, the plume 
of points turns off towards lower EWs, and the range 
and dispersion of EW exhibited by galaxies increases 
again. The sequence is then primarily occupied by bulge- 
dominated early-type spirals. 

We compute the changing dispersions and mean EW 
values along the sequence by collapsing Fig. la in four 
coarse bins of Mb chosen to correspond to the regimes 
defined by the two characteristic transitions (Table 1). 
Measurements of both Mb and EW are available for es- 
sentially all galaxies within the IIHUGS complete sub- 
sample, so the density of points in Fig. la is statistically 
representative of the Local Volume, down to our com- 
pleteness limit of —14.7. In all regimes, we find that the 
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TABLE 1 
EW Distribution Statistics 



EW(Ha+[NII]) 
no extinction correction 
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0.40 
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0.43 
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1.52 


0.16 


1.48 


0.17 


[-17.0,-14.7) 87 


1.50 


0.17 


1.49 


0.20 


[-14.7,-13.0) 102 


1.40 


0.33 


1.39 


0.35 



EWs are well characterized by log-normal distributions. 
For the two intermediate luminosity bins which include 
galaxies with —19 < Mb S —15, the means and disper- 
sions of the logarithmic EWs are essentially the same, 
with values of ~30A, and 1- and 3-a ranges (based on 
Gaussian fits to the distributions) from about 20-50A, 
and 10-lOOA respectively. Above the turn-off, there is a 
drop in the mean EW by a factor of ^2 to 15A, and an 
increase in the dispersion by a factor of 1.5. At the oppo- 
site end of the sequence, there is a factor of two increase 
in the dispersion in the lowest luminosity bin. There also 
appears to be a ~25% drop in the mean EW to 25A for 
these extreme dwarfs. We have checked that although 
the fractional errors in the EW do become larger with 
both decreasing EW and decreasing luminosity, changes 
in the average random uncertainty do not drive the vari- 
ation of the dispersion along the sequence. The average 
measurement error changes by several percent at most 
between luminosity bins, and this can cause an increase 
in the dispersion of only a few hundredths dex. 



4. DISCUSSION 

The plots in Fig. 1 provide a concise synthesis of the 
systematic dependencies of star formation activity with 
mass and morphology in the local universe. We com- 
ment on how these features reflect established empirical 
relationships, and re- frame open issues regarding star for- 
mation in sub-L* galaxies. 

First, we consider the galaxies that flank the transition 
at Mb ~ — 19 and Vmax ^ 120 km s^^ . Changes in the 
mean star formation history as a function of mass and 
morphological type has been well established for such 
systems (e.g., Kennicutt et al. 1994; Gavazzi & Scodeg- 
gio 1996), and the systematic drop in EW above the 
transition indicates that star formation activity is slow- 
ing down in galaxies above a certain characteristic mass. 
However, the abruptness of the transition to the "waist" 
of the sequence near Mb = —19 is somewhat surprising. 
This is not an artifact of our small 11 Mpc sample; the 
feature persists when we examine Ha EW datasets which 
probe bright galaxies over larger volumes (e.g. James et 
al. 2004, Moustakas & Kennicutt 2006). 

The narrowness of the sequence between — 19 < Mb ;$ 
— 15 (50 km s~^ < Vmax ^ 120 km s~^ ) places con- 
straints on the range of normal variation in star forma- 
tion activity as well as on the prevalence of starbursts. 
It implies both a high degree of temporal self-regulation 
within individual galaxies, and a small dispersion in evo- 
lutionary properties among the galaxies. We explore 
these issues in detail in Lee (2006) and Lee & Kenni- 
cutt (2008), but note here that based on the models of 
Kennicutt et al. (1994), the EW statistics imply that 
normal late- type spirals have b between 0.1 and 2 with 
an average value of ~ 0.5. That is, intermediate lumi- 
nosity disks commonly show factors of two to three fluc- 
tuations in their SFRs and are typically forming stars 
at half of their past average rate. In the currently fash- 
ionable downsizing paradigm, these are the galaxies with 
sufficiently low mass and high gas content that have not 
yet exhausted their interstellar gas supplies, and have 
not yet merged to form more massive systems. However 
such a scenario does not provide a ready explanation for 
why the star formation histories of such galaxies should 
be so homogeneous. 

Our results fit in well with recent work that has also 
aimed to characterize the star-formation properties of lo- 
cal galaxies. Brinchmann et al. (2004) used Sloan Digital 
Sky Survey (SDSS) spectra for ^ 10''' galaxies and exten- 
sive modeling to compute Ha-based SFRs. They show 
that b is nearly constant for 10^ < M^ < 10^°, but that 
the distribution broadens and shifts to lower values be- 
ginning at M, > lO^" Mq. Sahm et al. (2007) combined 
GALEX UV and SDSS optical photometry for -50,000 
galaxies and found that the UV-based specific SFRs also 
occupy a narrow sequence with mass and begin to decline 
on average in the same regime. Noeske et al. (2007) show 
that the population of galaxies with a narrow dispersion 
in SFRs at a given M, persists back to redshifts z ~ 1, 
and they refer to this feature as the "main sequence" of 
star-forming galaxies. Our new Ha narrowband imag- 
ing observations provide complementary confirmation of 
these results - they demonstrate that the trends also cor- 
relate with morphology and are already evident in the 
observable tracers of the physical quantities. 
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Work on understanding the physical origin of the tran- 
sition at Mb ^ — 19 and Vmax ~ 120 km s"^ (i.e., 
M, ^ 10^" 1^0 ) is on-going and generally involves inves- 
tigating gas inflow and/or energy feedback from super- 
novae and AGN into the ISM (e.g., Keres et al. 2005, 
Croton et al. 2006). Processes associated with bulge for- 
mation are also likely important, as we observe that the 
transition is accompanied by a shift from late- to early- 
type disks. We also note the intriguing coincidence of 
other disk properties which show transitions near Mb ~- 
19 and Vmax ~120 km s^^. The mass-metallicity relation 
undergoes a change in slope at approximately the same 
luminosity (Garnett 2002, Tremonti et al. 2004). Edge- 
on disk galaxies show distinct structures above and below 
Vcirc ~ 120 km s^^ , with more massive galaxies showing 
thin stellar disks and well-deflned thin dust lanes, while 
more slowly rotating disks exhibit thick stellar disks and 
no dust lanes (Dalcanton et al. 2004). This velocity has 
also been identified as the critical scale that subdivides 
systems that can retain the ejecta from supernova feed- 
back {Vcirc > 120 — 130 km s~^ ), as opposed to those 
that may lose some of their metals into the IGM (e.g., 
Martin 1999, Dekel & Woo 2003). 

Below Mb ^-15 and Vmax ^50 km s^^ the disper- 
sion increases again. With smaller representative sam- 
ples Kennicutt et al. (1994) and Hunter (1997) had al- 
ready noted that dwarf irregulars show a large scatter 
in their relative star formation properties as compared 
with spiral galaxies. In contrast to the upper mass tran- 
sition however, this feature has not been as well studied, 
because large surveys such as SDSS do not adequately 
sample the faint end of the galactic luminosity function, 
and previous studies of SFRs of dwarf galaxies lacked the 
completeness needed for statistical analysis. The lower 
transition is nearly as pronounced as the upper transi- 
tion, and it is interesting to ask whether it is indicative of 
another change in the dominant physical processes that 
regulate star formation. If so, then the transitions mark 
the boundaries of three regimes where distinct modes of 
star formation operate. 

A key question however is whether the increase in scat- 
ter at low masses reflects an increase in statistical fluctu- 
ations in massive star formation for galaxies with SFRs 
below ~0.01 Mq yr~^, rather than a fundamental physi- 
cal change in the degree of regulation in the global SFR. 
One possibility is that the increased scatter may arise 
from a breakdown in the reliability of Ha emission as a 
quantitative SFR measure when star formation is so low 
that it does not fully populate the upper end of the IMF, 
as suggested by Weidner & Kroupa (2006) and Pflamm- 
Altenberg et al. (2007) . However such effects are unlikely 
to account for the transition we observe. At Mb = —15 
the average SFR of the galaxies in our sample is approxi- 
mately 0.03 M© yr~^. A galaxy with the mean SFR and 
a Salpeter IMF would contain approximately 300 30 M© 
main sequence stars at any one time; if these stars formed 
independently, statistical fluctuations in their SFR of 
<10% would be expected, negligible in comparison to 
the observed dispersion in SFRs at any galaxy mass. If 
instead stochastic effects are dominated by the clustering 
of star formation, it is the \/N statistics of the number 
of clusters or HII regions that are relevant. However at 
the transition luminosity the average galaxy in our sam- 



ple possesses at least several HII regions, and often tens 
or more. Thus, this effect alone would only produce a 
dispersion of a few tens of percent, again too small to ac- 
count for the hundredfold range of SFRs observed. Still 
such stochasticity is important for the faintest galaxies 
in our sample (Mb > -13.5), with SFRs <0.001 M©, 
and we investigate this further in Tremonti et al. (2008). 

If instead the increased variability of the SFR at low 
galaxy masses is physical, then what mechanisms are re- 
sponsible? The answer to this question is still unclear 
although there has been much work on star formation in 
gas-rich dwarfs over the past three decades (e.g.. Hunter 
1997, Hunter & Elmegreen 2004, and references therein). 
The problem is particularly vexing because of the lack 
of an easily identifiable "second parameter" that cor- 
relates with the increased scatter; i.e., dwarf irregulars 
with very similar global properties (e.g., environment, 
gas content, structural properties) can have a wide range 
in normalized star formation activity (but see Papaderos 
et al. 1996, van Zee 2001). Our preliminary examination 
of EW with proximity to known neighbors, gas fraction 
and the concentration of Ha emission has also not un- 
covered any obvious trends. It is reasonable to speculate, 
however, that spiral structure plays a key role in regu- 
latory processes since the increase in scatter at the low 
end roughly tracks a shift in morphology from spiral to 
irregular types. 

The general absence of a relationship of activity in 
dwarf irregulars with environment outside rich groups 
and clusters (e.g., Lee et al. 1999, Hunter & Elmegreen 
2004) also suggests that internal processes may be re- 
sponsible for the increase in scatter. One possibility is 
that feedback from massive stars (e.g. stellar winds, su- 
pernova shocks) has a greater negative impact on the 
ISM of physically small, kpc-sized systems, with the re- 
sult that star formation is burstier in dwarf galaxies than 
in larger spirals. Gerola et al. (1980) originally explored 
such issues in their models of stochastic self-propagating 
star formation, which interestingly show a characteris- 
tic stabilizing transition in fiuctuations in the mean star 
formation rate for galaxies with Mb ^-15. This result, 
however, is dependent on the assumed "cell" size and the 
refractory period (the time during which negative feed- 
back inhibits the formation of new stars), which are free 
parameters in their toy models. The models of Mac Low 
& Ferrara (1999) incorporate more realistic prescriptions 
for the density distributions of the gas, stellar and dark 
matter in dwarf galaxies. Again, they find that below a 
rotational velocity of ^30 km/s, feedback disrupts and 
blows out the ISM from the galaxy. Most recently, Stin- 
son et al. (2007) have modeled the effects of supernova 
feedback in haloes with virial radius velocities from 10-30 
km s^^ using SPH + N-body simulations. These mod- 
els show there is an increase in episodic behavior with 
decreasing halo mass. 

This line of work offers a promising explanation for the 
low-mass broadening of the EW distribution, but a more 
careful comparison of model outputs and observations 
are needed. In particular, the combination of star forma- 
tion histories derived from observations of the resolved 
stellar populations of the nearest galaxies (e.g., from the 
ACS Nearby Galaxy Survey Treasury, ANGST, Dalcan- 
ton 2006) and the distribution of global properties from 
IIHUGS provides complementary constraints. Viable 
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models will simultaneously need to reproduce features 
of the temporally and spatially resolved star formation 
histories from ANGST, the trends in the relative SFR 
with mass as reported here, as well as the constraints 
on the starburst duty cycle which will be presented in 
upcoming papers. 
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